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Abstract 

In layered cuprates, the application of an in-plane magnetic field (H) changes 
the c-axis optical sum rule and superfluid density p s . For pure incoherent c- 
axis coupling, H has no effect on either quantities but it does if an additional 
coherent component is present. For the coherent contribution, different char- 
acteristic variations on H and on temperature result from the constant part 
(t±) of the hopping matrix element and from the part (t^) which has zero 
on the diagonal of the Brillouin zone. Only the constant part (t±) leads to a 

dependence on the direction of H as well as on its magnitude. 
PACS numbers: 74.20.-z, 74.25. Gz 
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A conventional c-axis optical sum ruletl has been observed in some high-T c cuprates such 
as optimally doped YBCOj! while in others, such as underdoped YBCO, it is violated. The 
sum rule is formulated in terms of the missing area at finite energy uj under the real part of 
the optical conductivity when the sample becomes superconducting. In the conventional case 
this missing spectral weight (AN) shows up as the superfluid density (p s ) at zero energy and 
AN/p s = 1. Underdoped YBCO exhibits a pseudo gap above the superconducting critical 
temperature (T c ). This is taken as an indication of non- Fermi liquid in-plane behavior and 
AN / p s ~ 1/2. Even for an in-plane Fermi liquid, deviations! from one result for AN/p s if 
the c-axis coupling is incoherent, i.e. does not conserve in-plane momentum. On the other 
hand, the conventional value of AN/p s observed in optimally doped YBCO is consistent 
with an in-plane Fermi liquid and coherent c-axis dynamics! whether or not the hopping 
amplitude depends on in-plane momentum. It is clear, therefore, that such measurements 
can yield important information on the dynamics of the charge carriers in the cuprates. 

Information on the nature of the interlayer coupling can also be obtained from other 
measurements. For example, it was recognized! that the temperature (T) variation of the 
c-axis penetration depth would mirror the observed in-plane linear T dependence, if the 
c-axis coupling is coherent with a constant matrix element t±. However, studies!'! of atomic 
overlapps between Cu02 planes through which c-axis tunneling occurs, have shown that the 
chemistry involved leads to a matrix element t(4>) which has (i-wave symmetry and is given 
by t^,cos 2 (20), where gives the direction of the momentum k in the 2-dimensional Cu02 
Brillouin zone. In general, there could be a contribution from both a constan t± and t(<f>). 
For the case of t(4>), the low T dependence of the resulting p s (T) for a d-wave superconductor 
is T 5 (See Ref.l). Such a dependence on T has been observed in a recent experimentB on 
Bi2212 indicating that t(<p) is the dominant coherent matrix element. While this experiment 
favors a T 5 , some small contribution from a constant t± and from some incoherent transfer, 
which gives a T 3 law,!'! cannot be ruled out. In fact, there is strong evidence from the c-axis 
quasiparticle tunneling work! on mesa junctions of Bi2212 that there is a significant t± part 
in the coherent matrix element .0 



In this paper we investigate the effect of an in-plane magnetic field H on the c-axis sum 
rule and on the corresponding superfluid density. The aim is to see if additional information 
on the nature of the c-axis coupling can be obtained from such experiments. We consider a 
general coherent c-axis coupling matrix element [ij_(k) = t ± + cos 2 (20)] as well as impurity 
mediated transfer from plane to plane. 

We begin with the Hamiltonian 7i = Hq + H c , where Hq describes a <i-wave su- 
perconductor in the plane and H c gives the c-axis coupling from plane to plane H c = 
So-,k,p^k-p C^ cr (k)C2a(p) + h.c. , where C^(k) creates an electron in plane i of momentum 
k and spin a, and tk-p is the hopping amplitude between plane 1 and 2. For coherent 
coupling, t k _ p = <5k- p t_i_(k). For impurity assisted hopping t k - P = Vk_ p , where Vk- P is the 
impurity potential, and in-plane momentum is not conserved once a configuration average 
over impurities is carried out. 

Consideration of c-axis electrodynamics leads directly to a relation for the superfluid 
density p s ; namely, p s = AN + K with the missing area AN = 8 Jq+ du [g™ c {oj) — cr{ c (u)] 
and K = —4ne 2 d [{H c ) s — (H c ) n ], where the cutoff (u c ) is to be applied at the energy of the 
interband transitions and o\l, s \u) is the real part of the c-axis optical conductivity in nor- 
mal(n) and superconducting(s) state, and e and d are the electron charge and the interlayer 
spacing, respectivly. The difference in c-axis kinetic energy between superconducting and 
normal state [(H c ) s — (H c ) n ] gives the correction to the conventional sum rule for which the 
kinetic energy difference vanishes. 

The effect of an in-plane magnetic field on the c-axis kinetic energy difference can be 
expressed in terms of the Green's function G(k, u) and the corresponding Gorkov function 
F(k,u). We assume that the field penetrates freely into the sample so that it is uniform 
between Cu02 planes. First, we consider the case of coherent c-axis coupling tj_(k). The 
c-axis kinetic energy is 

(# c )' = AT E ^( k ) 2 t G ( k + q> ^ G ( k ' £) - F(k + q, Q)F(k, Q)] 
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where G) = to + 7sgnu; with fermionic Mastubara frequency u, 7 is an effective impurity 
scattering rate modeled as a constant for simplicity, and q = (ed/2)z x H with z a unit 
vector along the c-axis. The Green's functions in Eq. (|1]) with k + q involve the quasiparticle 
energy -Ek+q — v/^k+q + ^ k with £k +q — £k + e q -^ Here, £k is energy spectrum measured 



with respect to the Fermi energy and e q = qvpcos((f) — 9 q ) with the Fermi velocity Vp and 
6 q the direction of q. Note that 9 q = 9 + 7r/2 can be interpreted as the direction of H (6) 
because of the symmetry in the problem. 

The difference (i/c) q — (H c ) s turns out to be, after some algebra, 
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The correction to Eq. (0) is O [(vFq/u c ) 2 (A /u c ) 2 }. Note that (vpq/uj c ) 2 is negligible and, 
furthermore, the quantity in the square bracket vanishes since T < w c . The same holds 
in the normal state i.e. (i7 c )q ~ (H c ) n — 0. Consequently, the kinetic energy difference 
between superconducting and normal state does not change with in-plane magnetic field. 
For incoherent coupling, the same conclusion applies but for simpler reason. In this case 

(77 c ) q ~ 4TiV(0) £ J ^^|^(0 fe ,0 p )| 2 | d£ k d£ p $ (e k + e q ,A k )$(e p ,A p ) , (3) 

where $(£, A) = (iu) + £ — A) / (cj 2 + £ 2 + A 2 ). By changing £ k + e q — > £ k , we see that 
(if c ) q = (H c ) s , and there is no change in kinetic energy. 

The superfluid density p s is p s = AN + if. In the presence of an in-plane field, the 
superfluid density (p s , q ) becomes AiV q + _K" q , where p S)q = p s + 5p s , AiV q = AN + <5AA^, 
and Kq = K + 5K. Since <5-fr = 0, we find Sp s = 6 AN. The conductivity sum rule (or the 
normalized missing spectral weight) is 

AiV q AN + 5 AN AN 5p 
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Ps,q Ps + 8p s p s P 

For pure coherent coupling, AN/ p s = 1, and for pure incoherent coupling, 5p s = 0; therefore, 
the sum rule is not changed by H. However, if both are present, the sum rule is changed 
because 5p s 7^ from the coherent contribution, and AN/ p s 7^ 1 from the incoherent part.i 



We next compute the change in superfluid density in the presence of H only for coherent 
coupling because H has no effect for incoherent coupling. For a constant t±, the superfluid 
density p Sjq is p Sjq = t 2 ± CT XL F(k + q, Co)F(k, Co), where C = 16ire 2 d. Assuming a 
cylindrical Fermi surface, we obtain, as T — > 0, 



27r e qV /4A2 + Cq 
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where = ^4A£ + e q ±e q and = */4A£ + e^A 2 + 7 2 ±7e q . At finite T with 7 = 0, 
we obtain 

k ^l-^k+q "^kJ 

where x{E) = tanh(E/2T)/E. For general coherent coupling, we simply replace t\ by 
£_i_(k) 2 within the sum over k. As H — > 0, p s , q (T) — > p s (T). Applying the nodal 
approximation^ to the general case for T < Ao, we obtain, at H = 0, for a pure case 
p s (T) = p s (0)-CN(0) [a x tl (T/A ) + a 3 t ± t^ (T/A ) 3 + a 5 t 2 (T/A ) 5 ] , where [a u a 3 , a 5 ) = 
[ln(2), 4.48^(3), 42.2£(5)]. The t\ contribution is linear in T as expected and the t 2 ^ goes 
like T 5 , a well known result.i The t^t^ cross term gives T 3 just as does the impurity medi- 
ated contribution.! In the above expression, p s (0) = CN(0) t\/2 + t±t ( f,/2 + 3^/16 . From 
here on we will use the notation p(* q with i — 1, 2, and 3 for t 2 , t±t^, and t 2 ^ contribution, 
respectively. 

Numerical results for p^q/P^(0) as a function of e q /T c with 9 q = and 7 = are given 
in Fig. (1). One should not confuse e q (= qvp) with e q = e q cos((p — 9 q ). The curve labeled (1) 
refers to the t\ contribution and is nearly a straight line implying a linear dependence on H 
(the magnitude of H). For t{4>) 2 = t 2 cos 4 (20), the effect of the field on pf\/pf\$) is now 
much smaller than for t\ as the curve (3) shows. This reduction results because the nodal 
quasiparticles are eliminated from the c-axis transport in this case by the factor cos 4 (20). 
The numerical data for the t 1 , case fit well an H 2 law at small H. The last curve in Fig. (1) 
which is labeled by (2) refers to the cross term 2tj_t^, cos 2 (20). It also fits an H 2 . Thus the 
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constant t± part in t±(k) is most effective in producing a change in p Sjq in the presence of 
H. 

It is instructive to derive an approximate, but analytic expression at T = 0, for the 
H dependence of p Sjq just discussed. We start from the temperature-independent part of 
Eq. (^). Applying a nodal approximation^ should be reasonable for the difference between 
p Sj q and p s (0). In this approximation, J2k ~^ ^node $ Jpdpd^i where J = {(2it) 2 vfVg] 1 
after an appropriate coordinate transformation. The integration over p is to be limited 
to the nodal region with cutoff po of order A . The quasiparticle energy = \f^~+ A 2 , 
then takes the form = Jp\ +pl = Pi where pi = pcos(-$) and P2 = psin(#). In J~, 
\/{tivfVg) = iV(0)/A o with N(0) the density of states at the FS. Near a node we can 
approximate e q ~ e q cos(0„ — 9 q ), where (fi n = ir/4, 37r/4, 57r/4, or 77r/4. Thus for a give 
value of 9 q , say 6* g = 0, e q ~ ±e q /\^2. 

For the case, choosing 9 q = 0, that is H along an anti-node, applying the nodal 
approximation we obtain p£ q /pW(0) = l-(7r/8 v / 2)(e g /Ao), where pW(0) = t 2 ± CN{0)/2. We 
choose A = 2.14 T c , then p^/p^O) = 1 - (e 9 /T c )/7.71. Comparison with our numerical 
result given in Fig. (1) shows that this simple expression works well up to e q = 0.5 T c . 
Similarly, for both t(0) 2 and t±t(4>) term, we obtain p^ q /p^(0) ~ 1 — (e g /A ) 2 /4, where 
i = 2 and 3. Consequently, in the presence of the in-plane magnetic field along the anti- 
node, the e q dependence of p Sjq is \e q \ for the t 2 ± term, and e 2 for both t±t(4>) and t(<p) 2 terms. 
Even if H is in the nodal direction, the e q dependence of p Siq does not change; however, the 
coefficient of |e g | (or e 2 ) becomes smaller because the field has no effect on the quasiparticle 
with (fi n ~ 9 ± 7r/2. 

While in Fig. (1) we took the field along the anti-nodal direction, in Fig. (2) we show 
additional results as a function of 9 in the range to ir for specific values of e q . The pattern 
obtained has four fold symmetry for a d-wave gap. Note that the underlying symmetry of 
the problem ensures that p s q is periodic in n/2 as is the quasiparticle conductivity.0 The 
thin solid line is for a constant t\ with e q = 0.1T C , 7 = 0, and T = while the thick solid 
line is for e q = 0.2T C (higher magnetic field). At all angle 9, p Sjq is reduced over its zero field 

6 



value. There is a minimum along the anti-nodal direction and a maximum along the nodal 
direction. In fact, for H along the node it can be shown by an application of the nodal 
approximation that p^q/p^(0) = 1 — (7r/16)e g /A . Note that the ratio of anti-nodal to 
nodal reduction in p Sjq is \[2 which is verified in our complete numerical calculations. This 
happens because the superfluid density depends on \e q \ for a constant hopping amplitude. 

Temperature and impurities will have an effect on the four fold symmetry pattern shown 
in Fig. (2). The first short dashed curve is for pJ>(T) with 7 = but now T = 0.001 T c 
while the second is for T = 0.01 T c . Both are to be compared with the T = result (thin 
solid curve). We see that in the first case T simply rounds off a little the maxima in the 
nodal directions while for a higher T = 0.01 T c , which is still much less than the magnetic 
energy e q — 0.1 T c , the temperature effects are more significant and much of the anisotropy 
is washed out. To see the anisotropy, which comes only from t\, one clearly needs to go 
to low temperature compared to the magnetic energy. Impurities also affect the anisotropy 
of Ps, q and results for 7 = 0.002 T c are shown as the long dashed curve in Fig. (2). These 
results are to be compared with the thick solid curve which applies for the same value of 
e q = 0.2 T c . Impurities reduces the c-axis superfluid density as we expect, and smears out 
the anisotropy. The final dot-dashed line in Fig. (2) is for the term t(4>) 2 with e q — 0.1 T c 
and 7 = 0. The superfluid reduction is small in comparison to the other curves and shows 
no anisotropy. 

More information on the effect of temperature, impurities and the magnetic field on the 
superfluid density for the t\ case is given in Fig. (3). All curves are for the field along the 
anti-node (except the solid curve in the lower inset, which is for the field along the node). 
The solid curve is for comparison and gives the T dependence of p^(T)/p^(0) with H = 0. 
The other curves are for different increasing values of e q /T c = 0.1 (dashed), 0.2 (dot-dashed), 
and 0.3 (long dashed curve). At low T, the power law in T is modified going to a T 3 while 
at high T it returns to T. It is seen that T always decreases the superfluid density as does 
e q . In all our finite T numerical results, we have approximated the T dependence of the gap 
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amplitude A(T) by Aotanh 1.74y (T/T c ) — 1 . In the top inset we show similar results 
but as functions of magnetic energy e q /T c for different values of 7/T c at T = 0. The solid 
curve is for 7/T c = 0, with which we compare other curves for 7/T c = 0.01 (dashed), 0.02 
(dot-dashed), and 0.03 (long dashed curve). Impurities, like T, have their greatest effect at 
e 9 = 0. 

We can obtain analytic results for the T dependence of p Sjq (T) = p Sjq + 8p SM (T). We 
start with the constant t\ case. Applying the nodal approximation already described, when 
the field is along an anti-node we obtain <5pW(T)/pW(0) = -12C(3) (T/A ) {T/e q ) 2 . For 
the t(<p) 2 term, we obtain 5 pi 3 q (T) / (0) ~ -9450C(7) (T/e q f (T/A ) 5 , which goes like 
T 7 . Similarly, we have, for t±t ((/>), 5 p^) (T) / pf ] (0) ~ -270C(5) (T/e q f (T/A ) 3 , which gives 
T 5 . For incoherent coupling, p Sjq (T) = p s {T) and its dependence remains T 3 as shown 
in Ref.i For coherent coupling, H along an anti-node changes the T dependence of the 
superfluid density as follows: p S)Q (T) ~ T 2 p s (T). If 9 ~ n , say 6 1 ~ 7r/4, then e q ~ for 
<fi n = 37r/4, 77r/4. This means that the in-plane field has no effect on quasiparticles with 
(fi n = 37r/4, 77r/4. Consequently, p Sjq (T) ~ p s (T) for H along a node. This is shown in the 
lower inset of Fig. (3). The solid curve is for the nodal direction while the dashed curve is 
along the anti-node. The anisotropy disappears for T > 0.1 T c . 

It is clear from our analysis that a contribution to p Sjq from a small t\ term will lead 
to anisotropy in p s q as a function of the orientation of H, and be linear in H for small 
field. In the pure case, its dependence is linear in T going over to T 3 at finite field along an 
anti-node with T < t q (magnetic energy). These characteristic variation are quite distinct 
from those predicted for a t((fi) 2 term or for the cross term t±t ((/>). In both these cases there 
is no anisotropy at T = and the changes with H are quadratic (ej?). In the pure case 
with zero field, the T variations are T 5 and T 3 , respectively, and go over to T 7 and T 5 at 
finite field for T < e q . These characteristic laws should help constraint further the relative 
magnitude of t±, and also of incoherent coupling. In this last instance, there is no effect 
of H and the T 3 remains at finite field and low T. 
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In conclusion, we have considered the effect of an in-plane magnetic field H on the c- 
axis sum rule and superfluid density(p S)q ). We have assumed that the in-plane dynamics 
can be described within a Fermi liquid theory. If the c-axis coupling is coherent, the sum 
rule remains conventional and unaffected by the application of H but p Sjq itself is reduced 
with increasing H. The reduction with H depends on the exact functional form of the c-axis 
coherent coupling and is different for a constant matrix element t± and t(<f>) = t$ cos 2 (20). In 
the first instance, to leading order, it varies like H (magnitude of H) while in the second like 
H 2 . If, however, the c-axis transport is incoherent, then both the sum rule and p Siq remains 
unchanged. When both coherent and incoherent interlayer transfer contribute, both the sum 
rule and p S;<i will change with H. The change in p sq comes only from the coherent part. 

For the case of coherent c-axis coupling with a constant t±, the change in p sq depends 
on the direction of H with respect to the nodes. For H along an anti-nodal direction p sq is 
minimum while it has a sharp maximum when H is along a nodal direction. This anisotropy 
in p sq is reduced when impurities are present or when the temperature is increased. For the 
case of t(<t>), the contribution to the transport coming from the nodes is eliminated and the 
anisotropy disappears. The incoherent coupling contribution also gives an isotropic response 
because in this case momemtum is not conserved and, correspondingly, the momentum shift 
caused by H is immaterial. 
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FIGURES 

FIG. 1. pi%/p^(0) as functions of e q /T c for T = 0, 9 = and 7 = 0. Note that the superscripts 

(1), (2) and (3) represent tj_, t±t<p and t^, respectively. 

FIG. 2. p^q(T)/pi^(0) as functions of 9. The dot-dashed line is for t\ with e q = 0.1T C , 7 = 0, 

and T = 0. The thin solid curve is for ^ with e q = 0.1T C , 7 = 0, and T = 0. The two dashed 

curves just below this curve are for t 2 ± with 7 = 0, and T = 0.001T C and T = 0.01T C , respectively. 

The thick solid curve is for the case of t\ with e q = 0.2T C , 7 = 0, and T = 0. The long-dashed 

curve is for e q = 0.2T C , 7 = 0.002T C , and T = 0. 

FIG. 3. pQ/pP (0) as functions of T/T c for different values of e q /T c with 9 = and 7 = (pure 

case). The solid curve is for e q /T c = (zero field). The dashed, dot-dashed, and long-dashed curve 

for € q /T c = 0.1, 0.2, and 0.3, respectively. In the upper inset, p^q/p^ (0) is given as functions 

of e q /T c for different values of 7/T c with 9 = and T = 0. The solid, dashed, dot-dashed, and 

long-dashed curve for 7/T c = 0, 0.01, 0.02, and 0.03, respectively. In the lower inset, phX/ p^P (0) 

as functions of T/T c for e q = 0.3 T c . The solid curve is for 9 ~ tt/4 and the long-dashed curve for 

9 = 0. 
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